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ABSTRACT 

We develop a parameter-free analytic model to include the effects of mergers into the theoretical modelling of the X-ray temperature 
function of galaxy clusters. We include this description into our model for the cluster population based on fluctuations of the gravi- 
tational potential, which avoids any reference to mass. Comparisons with a numerical simulation reveal that the theoretical model is 
in good agreement with the simulation results. We show that building the model on the dynamics of spherical rather than ellipsoidal 
collapse yields better results if emission-weighted temperatures are used, while ellipsoidal collapse yields good agreement between 
model and simulation for mass-weighted temperatures. Analysing two different samples of X-ray clusters, we quantify the influence 
of mergers and a conversion between different temperature definitions on the joint determination of CI^q and erg. If effects of mergers 
are included, temperature functions based on cluster masses and on the gravitational potential built on spherical collapse are in good 
agreement with other cosmological probes without any conversion of temperatures. 

Key words, cosmology: theory - methods: analytical - cosmology: cosmological parameters - galaxies: clusters: general - X-rays: 
galaxies: clusters 



1. Introduction 

Galaxy clusters are a potentially very powerful probe of non- 
linear cosmological structure formation since their abundance 
and its evolution depends sensitively on the matter density, 
the normalisation of density fluctuations and the dark energy. 
Conventionally, theoretical predictions of the cluster popula- 
tion parametrise clusters by mass. This is potentially problem- 
atic since mass is strictly not observable and an integral quan- 
tity which, for irregularly shaped bodies without well-defined 
boundary, is hard to define unambiguously. Calibration relations 
are needed between the mass and observable quantities such as 
X-ray temperature and luminosity, which are themselves prone 
to systematic and random uncertainties. 

We have proposed a different approach avoiding any refer- 
ence to mass (Angrick & Bartelmann 2009). The X-ray tem- 
perature function of the cluster population, i.e. their number- 
density distribution with X-ray temperature, can be theoreti- 
cally predicted based on the statistics of gravitational-potential 
fluctuations. This procedure has several advantages. First, it 
parametrises the cluster population directly by their temperature, 
which is a locally defined observable tightly related to the po- 
tential depth. Ambiguities caused by the integral definition of 
the mass are thus avoided. Second, calibration relations for the 
mass are circumvented, thus removing their scatter from the un- 
certainty of any inferences (see also Lau 201 1). Third, the grav- 
itational potential evolves much less than the matter density, ex- 
tending the range of validity of linear structure evolution. 

We have shown under which conditions this potential- 
based temperature function reproduces the theoretical predic- 
tions based on matter density and mass. Here, we address two 
subsequent questions. First, we compare the potential-based 
temperature function to a gas-dynamical, numerical simulation. 



While we find agreement at low redshift, there is increasing dis- 
agreement towards moderate and higher redshifts. This brings 
us to the development of an analytic model for the effect of clus- 
ter mergers on the X-ray temperature function, which leads to 
very good agreement of our theoretical predictions based on po- 
tential statistics with the numerical results. Our analytic model 
could be considered as providing an analytic complement to the 
numerical study by Randall et al. (2002). 

Second, we use the potential-based temperature function in- 
cluding the merger model to infer the cosmological parame- 
ters Qtno and erg from two different samples of galaxy clusters. 
The results are not conclusive yet, mainly because of tension 
between observationally infeiTed temperatures and theoretically 
motivated temperature definitions, but we find reasonable values 
for both parameters provided we use a definition of an X-ray 
temperature function that seems appropriate for the comparison 
with observational data. 

The paper is structured as follows: We first review briefly in 
Sect. 2 the derivation of the potential-based temperature func- 
tion, extending it to include ellipsoidal rather than spherical col- 
lapse. We compare it to a gas-dynamical numerical simulation 
in Sect. 3 and develop the analytic model for merger effects in 
Sect. 4. The inference of cosmological parameters is described 
in Sect. 5. Its results are discussed in Sect. 6, and we conclude 
with a summary in Sect. 7. 

2. The potential-based temperature function 

In the following, we use the X-ray temperature function intro- 
duced by Angrick & Bartelmann (2009) and an extension thereof 
based on the generalisation from spherical to ellipsoidal col- 
lapse. Both approaches avoid any reference to the globally de- 
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fined, strictly unobservable cluster mass, but are directly derived 
from the Gaussian statistics of cosmological potential fluctua- 
tions. 



2. 1. Original form: spherical collapse 

We briefly sketch the main idea and the basic steps in the deriva- 
tion of an X-ray temperature function for galaxy clusters that 
does not invoke cluster mass. It is based on the number density 
of minima of a homogeneous and isotropic Gaussian random 
field, discussed in great detail by Bardeen et al. (1986). For fur- 
ther detail on the derivation, we refer to Angrick & Bartelmann 
(2009). 

The differential number density of potential minima with 
depth between O and O + dO is 



n(<^) dO 



JaOc 



d(A<l))n(<D,A<l)) dO, 



(1) 



where the number density of relevant potential fluctuations 
n(<l), AcD) can be analytically expressed by 
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The quantity AO denotes the field's Laplacian and the lower in- 
tegration boundary is the critical Laplacian 



i5c(a) 



' ^mO - 



(5) 



where Hq is the Hubble constant, flmo is today's matter den- 
sity with respect to the critical density, a is the scale factor, and 
5^ is the critical linear density contrast of the spherical-collapse 
model. The spectral moments cr j of the redshift-dependent lin- 
ear potential power spectrum P<s,{k, z) are defined as 



(6) 



where W^{k) is the Fourier transform of a filter accounting for 
the shape of the gravitational potential of a homogeneous and 
isotropic overdensity. It is given by 



5 1^3 sin M - M (3 -H M^) cos mJ 



2m5 



(7) 



with u - kR and R - V-20/A0. The cut-off wave-number ^,„in 
introduces an effective sharp high-pass filter in A:-space and is 
chosen such that for a given combination (O, AO) the number 



density n(0, AO) is maximised. This filter step removes long- 
wave potential-fluctuation modes and thus ensures that clusters 
are not excluded that move along a potential gradient and thus 
do not have VO = at their centres. 

We have to relate the cluster's potential O to an observable 
quantity. This can be relatively easily done for its X-ray temper- 
ature in two steps: 

1 . The non-linear potential in the cluster centre is given by 

0„i = -2nGpi,6R\ (8) 

where G is Newton's constant and pb the background density 
of the Universe. It can be related to the linearly evolved po- 
tential using the spherical-collapse model for the evolution 
of the overdensity 6 and the radius R from an initial state 
at small scale factor to the time of collapse and final viri- 
alisation. The relation between the linear and the non-linear 
potentials, Oi and Oni, respectively, is then given by 



^ Rt. V'^ flc £>+(ai) 
Oi R^ C aiD+ia^ 



where 



C = 



^1/3 



1 + 



A,ta 



1 -Q. 



m,ta 



(9) 



(10) 



The subscripts 'i', 'ta', and 'c' refer to the initial, turn- 
around, and collapse times, respectively. D+ denotes the lin- 
ear growth factor of matter perturbations, and I is the over- 
density inside the sphere at the time of turn-around. Q,n and 
are the densities of matter and dark energy with respect 
to the critical density, respectively. 
2. The non-linear potential in the centre can be related to the 
local cluster temperature T using the virial theorem for 
Newtonian gravity, {Ey;-,^) = -\{E^ox), where (Ekin) = I^b^ 
and (/ipot) = -/^'WpOni are kinetic and potential energies, re- 
spectively, averaged over a sufficient number of particle or- 
bits. Here, yU « 0.59 (assuming that the intracluster gas has 
primordial composition with Helium abundance Y - 0.24), 
k^ is Boltzmann's constant, and is the proton mass. The 
central cluster temperature can then be expressed as 



1 



-jUmpOni. 



(11) 



The number density of galaxy clusters as a function of their X- 
ray temperature is finally given by 



Eq. (11) Eq. (9) 

n(T)dT^n(T % 0„i U Oi) 



dOi dOn 



dO„i dT 



dr. 



(12) 



2.2. Generalisation: ellipsoidal collapse 

The proposed temperature function is in good agreement with 
the mass function based on the Press-Schechter formalism (Press 
& Schechter 1974), but not with the Sheth-Tormen mass func- 
tion (Sheth & Tormen 1999, 2002) based on ellipsoidal-collapse 
dynamics, which better fits mass functions inferred from numer- 
ical simulations. Thus, we refined the relation between Oi and 
Oni by taking deviations from spherical dynamics into account 
as described in Angrick & Bartelmann (2010). In the following, 
we present the key ingredients and results of this model, which 
extends the original work by Bond & Myers (1996). 
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The dynamics of the dimensionless axes aj - Rj/Rp^^, j 
1, 2, 3 of a homogeneous elHpsoid, where /?pk the size of a spher- 
ical top-hat corresponding to a mass M 
scribed by the three coupled differential equations 



(47t/3)pb/?^^, are de- 
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with Cj = (1 + (5)/3 + bj/2 + /text,;- Here, E(a) is the expan- 
sion function of a Friedmann model, and the prime denotes 
the derivative with respect to a. The equations are coupled via 
6 = /(aia2ai) — 1, bj 
and external shear contributions with 



and /iextj- The latter two are the internal 
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and 



DM 



(linear approx.), 
(non-linear approx.), 



(15) 



respectively, where the subscript 'i' here and in the following de- 
notes the initial time. We use a combination of both shear mod- 
els, called the hybrid model, which describes the external shear 
by a linear approximation until the tum-around of an axis, and 
then switches smoothly to the non-linear approximation. 

The initial conditions for the evolution of the axes are derived 
from the Zel'dovich approximation and given by 



flj(fli) = fli[l — Aj{a{)\, 



dfl; 

dfl 



1 - 2^/fli). (16) 



The initial ellipticity e\ and prolaticity pi of the model are 
related to the initial overdensity 6i by 



3cr 



Pi =0, 



(17) 



where cr^ is the variance of the matter power spectrum. These 
values follow from the probability distribution of the eigenvalues 
of the Zel'dovich deformation tensor (Doroshkevich 1970). 

To stop the collapse, we use the following virialisation con- 
ditions for each axis, derived from the tensor virial theorem. 



fl'A 
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a^EHa) 



3^mO ^ 
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(18) 



The most important difference compared to the spherical col- 
lapse model is the circumstance that the parameters 5c and Ay be- 
come mass- or scale-dependent, respectively, and can be larger 
{6c for small masses) or smaller (Ay for large masses) by a factor 
of ~ 2 compared to the canonical values for the ACDM cosmol- 
ogy. We refer to Fig. 5 of Angrick & Bartelmann (2010) for their 
detailed dependence on mass and redshift. 

In the following, we will use the results of the ellipsoidal- 
collapse model and implement it in our formalism for the X-ray 
temperature function where results from the spherical collapse 
were used. 

We have to modify Eq. (1) since in the ellipsoidal-collapse 
case, the critical Laplacian AOc is now a function of the variable 



Aft, which one has to integrate over, through R = yl—2Q>\l h.<i> 
(see Sect. 2.1). It thus becomes 



n(Oi) d<l)] 



Jo 



d(A<E)) H^i, AO) e [A<D - AtDcCfi, AO)] dO,, 

(19) 

where is Heaviside's step function. Note that it is still a good 
approximation to smooth the density field with an isotropic top- 
hat of size R and not to introduce an anisotropic smoothing func- 
tion for the ellipsoidal-collapse model at the initial time, when 
we are well within the linear regime, and the deviation from 
sphericity is of order a few times 10"^. 

The potential in the centre of an ellipsoid is given by 



O = — 7TGpb&lfl2fl3^| 



pk 



Jo 



dr 



(20) 



fl? + T)(a? + T)(fl? + r) 



(Binney & Tremaine 1987, p. 57). For a sphere, fli = fl2 = 
fl3 = R/Rpk, and the integral can be solved analytically yield- 
ing 2Rp]^/R, hence the result for the sphere, Eq. (8), is repro- 
duced. We proceed exactly in the same way as in the previous 
section, calculating the ratio between linear and non-linear po- 
tential at the time of collapse. Again, quantities at a small ini- 
tial scale factor fl; are labelled with the index 'i' and quanti- 
ties at the time of collapse with 'c'. Using the approximations 
fli 1 ~ fli 2 ~ fli 3 ~ fli, and Ai * 1, we arrive at 



'Pnl _ flc 
Oi " 26i 



D+iflc) Jo 



dr 



a/(«Ic + ^)(«2,c + ^)(«3,c + ^) 



, (21) 



where we have also used the fact that for the virial overdensity, 
we can approximate (5v = Ay - 1 ~ Ay since Ay is of order 
100. All quantities that are necessary to evaluate Eq. (21) can 
be calculated using the ellipsoidal-collapse model by Angrick & 
Bartelmann (2010). Note that in the ellipsoidal case, the ratio of 
non-linear and linear potential becomes dependent on both Oi 
and AO. 

To infer an averaged linear potential for a given non-linear 
one, we marginalise over the dependence on AO weighted by 
n(Oi, AO) as follows. 



<Oi)A<I>(Oni) = 



/q" d(AO) Oi «(Oi, AO) e[AO - AOc(Oi, AO)] 

d(AO) «(4>i, AO) 0[AO - AOc(Oi, AO)] 

(22) 

where Oi - Oi(Oni, AO) via Eq. (21). 



3. Confronting theory with results from a simulation 

We now compare the analytic results for both X-ray tempera- 
ture functions, using the spherical- and the ellipsoidal-collapse 
dynamics, to a hydrodynamical simulation by Borgani et al. 
(2004) for a flat concordance ACDM model with matter den- 
sity Qn,() - 0.3, baryon density Obai-o = 0.04, Hubble con- 
stant Hq - 100 /i km s"' Mpc"' with h - 0.7, and normali- 
sation of the power spectrum erg = 0.8 in a box of side-length 
192 /i"' Mpc, starting at redshift Zstan - 46. The gas physics was 
implemented using gadget-2, a massively parallel A^-body/SPH 
tree code with fully adaptive time-resolution (Springel 2005). 
The density field was sampled with 480^ dark matter and an 
equal amount of gas particles with masses Mdm = 6.6 X 10'' Mo 
and Mgas = 9.9 x 10^ M©, respectively. During the time evolu- 
tion, the number of gas particles decreases due to their conver- 
sion into star particles, which have slightly lower mass than the 
gas particles. 
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The simulation includes radiative cooling processes follow- 
ing Katz et al. (1996), and a photo-ionising background expected 
from quasars which ionise the Universe at z ^ 6. Star forma- 
tion is modelled using the hybrid multiphase model for the inter- 
stellar medium by Springel & Hemquist (2003). The simulation 
code also includes a method to follow the production of metals. 
However, the effects of metals on the cooling function are not 
taken into account. This only affects the analysis of simulated 
galaxy clusters with temperatures T < I keV. For the evaluation 
of the X-ray temperature function, cluster catalogues for the five 
redshifts z = 0, 0.2, 0.5, 0.7, 1 are available. 

3. 1 . Problems in the analysis 

In Sect. 2, we established a relation between the depth of a clus- 
ter's potential minimum and its temperature in the centre. In the 
derivation, however, we neglected additional difficulties arising 
from baryonic physics, such as additional cooling of the gas and 
feedback processes, e.g. from supernovae or AGN. In combina- 
tion, these effects lead to cool cluster cores. These are in direct 
contrast to the naive expectation from Eq. (11) that the temper- 
ature in the centre should have a distinct maximum. Moreover, 
the temperature inferred from the measured cluster spectrum dif- 
fers from the temperature to be inserted into the virial theorem 
mainly for the following reasons. 

1. A detector counts the number of photons emitted by the 
source as a function of their energy. Hence, the more pho- 
tons are counted in an energy bin, the larger is its weight 
in the determination of a cluster's temperature. Parts of the 
cluster having different temperatures are thus weighted by 
their emissivity. 

2. The detector response is not the same for all energy bands. 
Although the numbers of photons collected in two energy 
bins might be equal, their weight will be different. The de- 
tector's sensitivity needs to be accounted for 

These two aspects are reflected in the different temperature def- 
initions that are used in numerical simulations. There, the tem- 
perature of a cluster is defined as 

f wr,as dV 

T = , (23) 

/ WAV 

where Tgas is the temperature of a gas particle and W is a weight 
function (Mazzotta et al. 2004). The integral covers a specified 
volume, e.g. the sphere defined by the virial radius. Depending 
on the choice of the weight function, mainly three different tem- 
peratures are used in the literature: 

1. the mass-weighted temperature T^^ with W - Wgas, where 
Hgas is the number density of the gas, 

2. the emission-weighted temperature Tew with W - 
MTgna) Hgas' whcrc A(rgas) oc ^JT^ IS the cooling function, 

3. the spectroscopic-like temperature with W - nl.^^/Tll^. 

Tjnw is easy to calculate and the temperature to be used in the 
context of the virial theorem, Eq. (11). However, it differs from 
the temperature inferred from observations. Consequently, Tew 
was used to relate temperatures from simulations with spec- 
troscopically derived temperatures. But since disagreements re- 
mained, especially due to the inhomogeneous sensitivity of a de- 
tector across its energy bands, was introduced to match the 
spectroscopic temperatures of clusters in Chandra and XMM- 
Newton surveys better. 



For z = 0, the catalogue includes three differently defined 
temperatures averaged over various radii. The same is true for 
the higher redshifts, besides that is missing. 

Also based on the simulation by Borgani et al. (2004), Rasia 
et al. (2005) provide a fitting formula relating the spectroscopic- 
like to the emission-weighted temperature for clusters with 
Tew > 2 keV at z = 0, 

r^i = (0.70 + 0.01)re„ + (O.29 ± 0.05). (24) 

We will use that relation in Sect. 6 later 

In addition to the aforementioned differences in the temper- 
ature definitions, hydrodynamical simulations suffer from the 
overcooling problem. The cooling function used in simulations 
is too efficient in the centres of clusters compared to observa- 
tions so that more gas can fall into the inner regions than what 
should be actually allowed. On its way, the gas is adiabatically 
compressed and heats up. Simulated cluster centres tend to have 
much higher temperatures than those observed (see Borgani & 
Kravtsov 2009, especially their Fig. 9). 

A more careful analysis of how to compare properly temper- 
atures that are inferred from both simulations and observations 
to those used in our theoretical model remains to be done. We 
postpone it, since first, it requires a detailed review of each indi- 
vidual step in the X-ray analysis of observed clusters, and sec- 
ond, the cluster temperature definitions in the theoretical work 
might also have to be altered, yielding the possibility to compare 
them to temperatures inferred from numerical simulations with- 
out suffering too much from the overcooling problem. A thor- 
ough study of these points will be provided in a forthcoming 
paper. 

3.2. Comparing different temperature functions 

Comparing the different temperature definitions averaged inside 
the virial radius /?vii with those averaged within i?25oo, where 
^2500 = (0.12 - 0.33) /?vii, i e- only within the inner part of 
the cluster, where the overcooling problem in the simulations 
occurs, yields interesting results. While the emission-weighted 
temperatures inferred by averaging over the two radii are essen- 
tially the same, both the mass-weighted and the spectroscopic- 
like temperatures differ within /?vir and ^2500- At z = 0, 
7'raw(^250o) is -45% larger than T^s^iR^i,), and T^i(R250o) is 
~20% larger than T^iiR^ij), showing that is most severely 
biased by overcooling. 

In Fig. 1, we compare the differential number density of clus- 
ters inferred from the simulation for z = using all three dif- 
ferent temperature definitions, where the theoretical prediction 
includes either spherical- or ellipsoidal-collapse dynamics. The 
temperature function predicted based on ellipsoidal collapse is in 
better agreement with the simulation if the mass-weighted tem- 
perature within the virial radius /?vir is used (upper panel), while 
incorporating spherical collapse leads to a theoretical tempera- 
ture function that agrees well with the simulation if emission- 
weighted temperatures are used instead. 

The reason for this discrepancy is that T^^ is weighted by nl.^^ 
and therefore accentuates the innermost cluster cores which ex- 
perience overcooling, whereas r^w is weighted by «gas, whence 
outer parts become more important. The overcooling problem 
is thus most severe for the temperature function based on Tew. 
In this case, the number density of clusters based on the core 
temperatures is only accidentally well modelled with spherical- 
collapse dynamics. The weight for the temperature function 
based on T^i causes it to fall between the two predictions. 
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Fig. 1. Comparison of the X-ray temperature function inferred 
from the simulation at z = using different temperature def- 
initions with the theoretically predicted temperature using ei- 
ther the spherical- (Osph) or the ellipsoidal-collapse model (^eu). 
Upper panel: Temperatures evaluated inside Ry„. Lower panel: 
Temperatures evaluated inside 7?2500- 



Interestingly, when considering the temperature functions in- 
ferred within /?2500 (lower panel), i.e. restricting the point of view 
to the very inner part of the halo per definition, the three func- 
tions approach each other and are more or less well described 
by the theoretical description based on the spherical-collapse 
model, where again the temperature function based on Tew fits 
the prediction best. This behaviour agrees with our earlier re- 
mark that the inner parts of clusters that experience overcooling 
are well described by spherical-collapse dynamics. 

Since T^i^ is the temperature that should be used with the 
virial theorem Eq. (1 1), we concentrate on the temperature func- 
tion based on T^w within 7?vir from now on. Figure 2 compares 
the results of the simulation based on T^^ with the potential- 
based temperature function including ellipsoidal collapse. One 
can clearly see that the higher the redshift is, the more the re- 
sults of the simulation disagree from the theoretical prediction. 
The simulation data systematically fall above the prediction for 
z > 0.5. In the following section, we hypothesise that this is 
due to cluster mergers, which are much more numerous in the 
past, and develop a model for taking them into account in the 
construction of the cluster temperature function. 




T [keV] 




Fig. 2. Comparison of the X-ray temperature function from the 
simulation based on T^w within /?vir with the theoretical predic- 
tion including the ellipsoidal-collapse model. 



4. Modelling merger effects 

Randall et al. (2002) have shown with numerical simulations that 
mergers do have a strong impact on the X-ray temperature and 
luminosity functions of haloes and thus, cosmological parame- 
ters infeiTed from them without including this additional effect 
are biased. Since clusters that are undergoing mergers are shifted 
from lower to higher temperatures and due to the exponential 
cut-off" at the high-temperature end, especially this part of the 
temperature function is enhanced. Consequently, Randall et al. 
(2002) find that the infeiTed erg is biased towards higher values, 
whereas Q,no is biased towards lower values. Both parameters 
change by ~ 15-20%. 

We choose a diff'erent approach here, trying to incorpo- 
rate the essential physical eff'ect of mergers in a simple ana- 
lytic model based on the extended Press-Schechter formalism 
by Lacey & Cole (1993). They found that the conditional proba- 
bility for a halo with mass M\ at time t\ to have mass M2 + AM2 
at f2 is given by 



^(52,W2|5 1,(^1) d52 = 



1 



Si 



V27T 

X exp 



3/2 



<ji)2{(JL>\- CJ2) 



{U)2S\ - Ci)\S2f 



2S,S2{S,-S2) 



dS2, 



(25) 



where S 1 and S 2 are the variances of the linear matter power 
spectrum filtered on scales corresponding to masses Mi and M2, 
respectively. The w,, / = 1,2, denote the scaled critical linear 
density contrasts from the spherical-collapse model at times f, 
defined as w,- = 6c,i(ti)/D+{tj). Note that > ^2 and oji > (±)2 
if defined in this way. Changing variables to M = M\, AM = 



5 



C. Angrick and M. Bartelmann: Cluster mergers and cosmological parameter constraints 
M2 - Ml and z = zicoi), Az = z(coi) - zicoi) yields 



temperature 

T 



p(M, AM, z, Az) d(AM) = p[S (M + AM), a)(z)\S (M), aj(z - Az)] 

dS 



X 



d(AM) 



(M + AM) 



d(AM). 



(26) 



Assume now that the temperature increase AT(M, AM) due 
to a merger of a mass M with another mass AM < M originates 
from the kinetic energy of the gas transported with the infalling 
clump, which is completely transformed to thermal energy. The 
gain of energy is therefore 



AE = ^/bAM(Av)2 i ^Nk^AT, 



(27) 



where f^AM is the baryon fraction of the lower-mass halo, Av 
is the relative velocity of the components, and = fbM/ifinij,) 
is the total number of gas particles in the halo of mass M. Note 
that in this ideal case, the factor /b cancels exactly. To guess 
Av, assume that the larger component is at rest, while the halo 
of mass AM approaches from infinity. In this case, the velocity 
can be easily calculated by equating potential and kinetic energy. 
The potential at the surface of a homogeneous ellipsoid is given 
by 



3MG 



4 R 



pk 



E\{a\,a2, 03) - ^ x^j E2{a\,a2, ay,,aj) 



with 



f" dr 

Ml, 02,03,, an - I — z 

' Jo (fl +T)nli(«^+T)i/2 



(28) 

(29) 
(30) 



(Binney & Tremaine 1987, p. 57), where we have used that the 
mass of the ellipsoid is given by 



47T 

M = — pa ifl2a3-'?pk 



(31) 



(cf. Sect. 2.2), and x denotes a position on its surface. A prop- 
erly averaged potential (O), where (■) denotes an average over 
all directions, can be calculated by introducing ellipsoidal coor- 
dinates xi — ai sin 6 cos (p, X2 — 02 sin 6 sin (p, X3 = 03 cos and 
averaging over the two angles 6 and (p. This finally yields 



(0) = - 



3MG 



4 R 



pk 



£■1(01,02,03) - T ^a^jE2{ai, 02,03, aj) 



Since Av = V=2<0), Eq. (27) yields 

2 pm^ AM<<1)) 



M 



(32) 



(33) 



where (O) is given by Eq. (32). We assume that the time scale 
for the temperature increase is set by the sound-crossing time 



R 



with 




(34) 



(Randall et al. 2002), where c^ is the sound speed. To infer 
a proper radius R, assume that the halo's mass is given by 



temperature boost 


due to mergers 


%-„_i(M) + AT{M, AM) 






virialised object 






T„,_2{M + AM) 


virialised object 













time t 



Fig. 3. Illustration of a cluster's temperature curve due to a 
merger as assumed in our model. 

M = (47t/3)pbAv^'^, where Ay - 0^/(010203) in the ellipsoidal- 
collapse model, so that 



R(M) = 



(3M01O203 
47tpcr0f^i 



■mO 



1/3 



(35) 



where p^o is the critical density of the Universe today. 

As illustrated in Fig. 3, we thus model the temperature boost 
in an idealised, abrupt way: The halo of mass M has a tempera- 
ture rvir,i(M) before the merger, which increases instantaneously 
to Tyii iiM) + AT{M, AM) for a period fjc and then drops instan- 
taneously to Ty„2{M -H AM), assuming that a virialised halo of 
mass M + AM has finally formed. To assign a virial temperature 
Tvii- to a halo of mass M, we construct a relation that is based 
on the combination of Eqs. (11), (20), and (31), thus assuming a 
virialised homogeneous ellipsoid. This yields 



k-RT^i,(M) = pm^ 



128 



dr 



» Uli[al(M) + T] 



1/2- 



(36) 

Starting from the number density of virialised galaxy clusters 
rinriT), which can be calculated as explained in Sect. 2.1, we 
calculate two correction terms. The first is the number density of 
clusters that reach a temperature T only due to mergers. 



n+(r)s I d(AM) I dMn[T^ir(M)]p{M,AM,z,Az) 
Jo Jo 

xSd[T- T^iAM) - AT(M, AM)] e(M - AM), 



(37) 



where we ensure via Dirac's delta function 5d and Heaviside's 
step function 9 that only combinations of AM and M contribute 
to the integi-al for which rvir(M)+ Ar(M, AM) = T and AM < M 
are fulfilled. The redshift interval Az is set by the sound-crossing 
time, Eq. (34), since in our simple model, one should be able to 
see all mergers at redshift z that have happened in the redshift 
interval Az before. The relation to fsc is given by 



c+Az 



dz' 



E{z')(l+z') 



(38) 



Two assumptions where implicitly made during the deriva- 
tion: First, the number density n{T) is assumed not to change sig- 
nificantly during Az, and second, the increase by AM is only due 
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Fig. 4. Comparison of the result from the simulation with the 
theoretical prediction both including (solid line) and excluding 
(dashed line) the effect of mergers on the temperature function. 
We compare with the temperature function derived from the sim- 
ulation that is based on within 7?vh (red data points). 

to a single merger event, ignoring multiple simultaneous merger 
events and smooth accretion. It turns out, however, that Az is 
short enough for these assumptions not to result in a significant 
error contribution. 

The second correction term arises due to clusters that would 
have a temperature T if they were virialised, but have a temper- 
ature higher than T due to mergers. 



n-(T) = 



M{T) 



d(AM) n(T)p[M(T), AM, z, Az]. 



(39) 



Here, we have to invert Eq. (36) numerically to assign a mass 
M to the temperature T. The total number density of clusters is 
then given by 



«ges(r) = n^iriT) + n+{T) - n_{T). 



(40) 



In Fig. 4, we compare the prediction of our X-ray tempera- 
ture function both with and without the enhancement by merg- 
ers with the results from the simulation by Borgani et al. (2004). 
The shape of the temperature function influenced by mergers is 
in qualitative agreement with Randall et al. (2002). The relative 
change of the number density increases with temperature and 
redshift. Using our simple model for the effects of mergers, the 
temperature function based on the mass-weighted temperature 
within 7?vii from the simulation is now in good agreement with 
the theoretical prediction. 



temperature function both with and without the effect of mergers 
to quantify their influence on the final outcome. Additionally, 
we compare the results of our approach using the statistics of 
minima in the cosmic gravitational potential to the traditional 
method invoking mass functions and an empirical M-T relation 
to see if any of the two gives tighter constraints. 

5. 1 . The samples 

The first flux-limited sample by Ikebe et al. (2002) consists of 
61 clusters and is based on ASCA and ROSAT data with a me- 
dian redshift of z = 0.046 in the temperature range 1 .4 keV < 
T < 10.55 keV. It covers 8.14 steradians, and the flux limit is 
1.99 X 10"" ergs s"' cm"^ in the 0.1-2.4 keV band. The maxi- 
mal search volume y^ax for each cluster is calculated and listed 
for an open model with Q^o = 0.2 and Qao = 0.0 and for a 
flat model with Q^o = 0.2 and Q^q - 0.8. Although Vi^^x itself 
depends on the cosmological parameters, it changes only very 
slightly with them. Neglecting this eff'ect in the following analy- 
sis therefore does not induce a significant eiTor. 

The second sample encompasses two subsamples by 
Vikhlinin et al. (2009a), one at high and one at low redshift, 
based on ROSAT PSPC All-Sky (RASS) and 400 deg^ data. 
The low-redshift sample consists of 49 clusters with flux / > 
1 .3 X 10"' ' erg s"' cm"^ in the 0.5-2 keV band from several sam- 
ples of RASS with a total area of 8.14 steradians. The redshift 
coverage is 0.025 < z < 0.25 with (z) ~ 0.05, and temperatures 
are in the range 2.61 keV <T < U.ll keV. 

The high-redshift sample consists of 36 clusters from the 
ROSAT 400 deg2 survey (Burenin et al. 2007) in the redshift 
range 0.35 < z < 0.9 with (z) ~ 0.5 and a redshift-dependent 
flux limit in the 0.5-2 keV band. For z > 0.473, the limit- 
ing flux is 1.4 X lO"'-' erg s"' cm"^, whereas for z < 0.473, 
the flux limit coiTesponds to a minimal X-ray luminosity of 
Lx,min = 4.8 X 10''^(1 + z)'-^ erg s"'. The temperatures of the 
clusters are in the range 2.13 keV < T < 1 1.08 keV. 

For both subsamples, the effective differential search volume 
dy/dz as a function of mass M and cosmological parameters 
QmO, f^AO and h - 0.72 for both subsamples was made avail- 
able in electronic form on a grid by A. Vikhlinin. To convert it 
to a function of temperature, we used the best-fit values of the 
mass-temperature relation of Vikhlinin et al. (2009a), 



M500 = M( 



\5keV/ 



(41) 



where Mq = (3.02 + 0.11) x lO''* /i"' Mq, a = 1.53 + 0.08, and 
E{z) is again the expansion function of a Friedmann model. 

5.2. The fitting procedure 

Since the errors on the cluster number counts are Poissonian, we 
use the C statistic of Cash (1979) for unbinned data to find the 
best-fit values for Q^o and erg, assuming a spatially flat universe, 
hence Qao = 1 - Qmo- The C statistic is defined as 



C = 2 



A^- y Inn, 



(42) 



5. Inferring cosmological parameters 

In this section, we want to infer the cosmological parameters 
Qmo and (Tg from two different samples by Ikebe et al. (2002) and 
Vikhlinin et al. (2009a) using our theoretical model for the X-ray 



where is the total number of objects expected from the sample 
assuming a theoretical model and is the theoretically expected 
differential number density of the /-th cluster in the sample with 
temperature T, and redshift z,. The sum extends over all sample 
members. 
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Although the potential-based temperature function using 
spherical collapse seems to be in good agreement with the tem- 
perature function from the simulation based on Tew by Borgani 
et al. (2004) only because of the overcooling problem, we also 
include it in our further analysis and either apply a temperature 
conversion according to Eq. (24) or identify the measured tem- 
perature directly with the one used in the theoretical framework. 
In these cases, we simply set ai - 02 = 03 = R/Rpk in our 
merger model. 

Thus, we fit in total eight different theoretical models to the 
two subsamples of Vikhlinin et al. (2009a): 

1. mass function by Tinker et al. (2008), 

2. the same including merger effects, 

3. temperature function based on the gravitational potential 
including spherical-collapse dynamics without temperature 
conversion, 

4. the same including mergers eff'ects, 

5. temperature function based on the gravitational potential in- 
cluding spherical-collapse dynamics with temperature con- 
version, 

6. the same including merger effects. 

7. temperature function based on the gravitational potential in- 
cluding ellipsoidal-collapse dynamics without temperature 
conversion, 

8. the same including merger eff'ects. 

In the first two cases, we assume a mass-temperature rela- 
tion according to Eq. (41) whenever we have to relate a mass to 
a temperature or vice versa, thus also when applying our analyt- 
ical merger model. To properly take the scatter into account, we 
convolve with a log-normal distribution. 



T=5 keVx[M E(z)/Mo]^ 



PMiT\M)dT^ 



1 



27TO-in7-r 



exp 



[T - Tq(M )] 



2^ 



dr (43) 



where To(M) is given by Eq. (41) with Mq = 3.02 x 10'^ h'^ Mq 
and a = 1.53. The standard deviation is mass-dependent and 
given by 




25 



Mo=3.02x10''h ' M^^n, a=1.53 
a: ±0.08 
O|m-=0.03+0.041ln[M E(z)]-ln MqI 



M E(z) I 



=(T„-a)/p 




T31 [keV] 

Fig. 5. Uncertainties in the mass-temperature relation {upper 
panel) and the Ts\-Tg„ relation (lower panel), modelled via vary- 
ing standard deviations of a log-normal (crj,, 7-) and a normal dis- 
tribution (cTj), respectively. 



cr,„7- = 0.03 + 0.04 \ln[ME(z)] - In MqI 



(44) 



due to the uncertainty of Mo and a; see the upper panel of Fig. 5. 
In cases 5 and 6, we convolve with a normal distribution of the 
form 



PTAT\Tsi)dT 



1 



V27t 



exp 



[T - rew(r ,i)] 

2cr2 



2\ 



dT, (45) 



where T^^ = (r^, - 0.29)/0.7 (Eq. 24). To take the scatter in 
the relation into account, the standard deviation is set to ctt - 
0.1 V^d^; see the lower panel of Fig. 5. We also use Eq. (45) to 
account for an error contribution in the remaining cases 3, 4 and 
7, 8 by simply setting Tsi = Tew. 

Taking these uncertainties into account, the expected number 
of objects of each subsample is given by 



dT^-^^(T,z) 
dz 



dxn(x)px(T\x), 
(46) 

where x can be either mass M or spectroscopic-like temperature 
Tsi depending on the theoretical model used to fit the data (mass 
function or temperature function). The integral boundaries de- 
pend on the subsample and are given in Sect. 5.1 for z and T. 
The integration over x has to be done over the whole valid range 



of p;t. Finally, the expected differential number density of the 
i-th cluster is simply given by the convolution 



n 



/,low|high ~ 



dVlowlhigh 

"d^ 



,Zi) J dxn 



(x)pATi\x). (47) 



To jointly fit both the low and the high-redshift cluster samples 
of Vikhlinin et al. (2009a), we have to add the two contributions, 
finding 



Cv = 2 



Mow - In n/,iow + Nugh - ^ In n 



/.high 



(48) 



For the sample by Ikebe et al. (2002), we proceed analo- 
gously, but the situation is much easier since we only deal with 
one single sample that covers only a small redshift interval. The 
latter implies that we do not introduce a significant error if we ig- 
nore the redshift evolution of the mass or temperature function, 
respectively, in the analysis. Instead, we compute the theoreti- 
cal functions at the sample's median redshift of z = 0.046 in 
the same way as Ikebe et al. (2002) did so that we only have 
to integrate over the temperature when calculating the total ex- 
pected number of objects from the sample. Hence, the C statistic 
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is given by 



C, = 2 



(49) 



with 



N = j^'dT V^,,(T) j dxn(x)p,(T\x), (50) 

«/ = Vn,ax(7'/) J dxn(x)pATi\x). (51) 

The conditional probability p^{T\x) is again given by Eqs. (43) 
or (45), respectively, thus assuming the same errors on the rela- 
tions as for the Vikhlinin et al. (2009a) data. To better compare 
with the results by Ikebe et al. (2002), we shall also use the clas- 
sical Press-Schechter mass function instead of the one by Tinker 
et al. (2008) and relate mass and temperature via Eq. (36) with 
ai - a2 - 03 = R/Rpi^. 

We search for minima of the C statistic as a function of the 
two cosmological parameters Q^o and erg, which enter both via 
n(x) and the volume factors dV/dz and Vmax- Only because the 
latter is very insensitive to changes in these two parameters, 
its dependence on Q^o and erg can be neglected. Cash (1979) 
showed that one can create confidence intervals for the C statis- 
tic in the same way as it can be done for ax^ fit using properties 
of the;^'^ distribution. Following Lampton et al. (1976), intervals 
with confidence y are implicitly given solving 



y 



-f 

Jo 



(52) 



for f, where / is the density of the distribution with p de- 
grees of freedom determined by the number of parameters. For 
95% confidence and p - 2, it follows that t - 5.991. Using the 
minimum of the C statistic, Cmi„, we can simply calculate the 
95% confidence contours by searching for points in the parame- 
ter space for which C - Cj^in + 5.991. 

6. Results 

In Fig. 6, we present the 95% confidence contours for both the 
Ikebe et al. (2002) and the Vikhlinin et al. (2009a) samples. 
Comparing the upper and the lower panels, one can see that the 
results from both data sets are compatible with each other, al- 
though pronounced differences exist. However, for the latter, the 
confidence contours are smaller due to the additional informa- 
tion on the redshift evolution of the temperature function. 

The mass-based temperature functions (red contours) and the 
potential-based temperature functions using spherical collapse 
without temperature conversion (blue contours) give similar and 
compatible results, i.e. the confidence contours are in agree- 
ment with each other and have approximately the same size. 
When redshift evolution information is added (lower panel), it 
seems that the direction of degeneracy is slightly changed for 
the potential-based temperature functions so that erg can be con- 
strained more tightly compared to the upper panel. For the mass- 
based temperature function, there is no such effect. Overall, the 
potential-based functions yield slightly smaller best-fit values 
for cTg compared to their mass-based counterparts. 

If merger eff'ects are taken into account, erg is significantly 
lowered while Qmo is increased for all temperature functions. 
Consequently, the confidence contours of the results including 
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Fig. 6. 95% confidence contours for the various theoretical mod- 
els presented in Sect. 5.2. The results for the potential-based 
temperature function assuming spherical-collapse dynamics are 
labelled with Osph, those incorporating ellipsoidal collapse with 
Oeii, respectively. While the solid curves represent the results 
of temperature functions that take into account merger effects, 
the dashed curves are for those which do not account for merg- 
ing. The crosses mark the respective best-fit values. Parameter 
constraints from a joint analysis of WMAP5, BAO, and SN 
data inferred from Komatsu et al. (2009) are indicated by the 
black solid contour. Upper panel: Ikebe et al. (2002) data. Lower 
panel: Vikhlinin et al. (2009a) data. 



and excluding mergers do not overlap. This result is in good 
agreement with the work by Randall et al. (2002), who also 
found using numerical simulations that mergers do have a drastic 
impact on the results for those two parameters. 

First converting the measured temperature according to 
Eq. (24) and then comparing it to the theoretical potential-based 
temperature function incorporating spherical collapse (green 
contours) has a similar effect as mergers: erg is decreased while 
is increased. Interestingly, using the potential-based temper- 
ature function including merger effects but without temperature 
conversion (blue solid contours) and using the potential-based 
temperature function with temperature conversion but without 
merger effects (green dashed contours) give almost identical (as 
for Vikhlinin et al. 2009a) or at least similar (for Ikebe et al. 
2002) confidence contours. Hence, the effects of temperature 
conversion and cluster mergers turn out to be highly degener- 
ate. Additionally, merger boosts increase the uncertainty on Q^o 
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since the confidence contours become more elongated in this di- 
rection. 

The confidence contours inferred from the potential-based 
temperature functions built on ellipsoidal-collapse dynamics 
(magenta contours) are shifted towards higher Q,no compared to 
those built on spherical collapse, both including and excluding 
mergers. This can easily be understood considering Fig. 1 again. 
There, the temperature function based on ellipsoidal collapse lies 
always below the spherical-collapse case for the same cosmolog- 
ical parameters. The reasons are a significantly enlarged and a 
smaller Ay that leads to a smaller ratio of the non-linearly to the 
linearly evolved potential. The result of the former is a decrease 
of the temperature function's amplitude while the latter shifts the 
curve towards smaller temperatures. Hence, to arrive at a similar 
fit for a given data set, a higher flmo is needed. Additionally, the 
contours based on ellipsoidal collapse are more extended espe- 
cially in the Omo-direction compared to both the spherical case 
and the mass-based temperature function. This implies that the 
potential-based temperature function built upon ellipsoidal col- 
lapse is the least sensitive to changes in 0.^0- 

The contours are not in agreement with a joint analysis of 
the five-year data release of the Wilkinson Microwave Anisotropy 
Probe (WMAP5), and data from baryonic acoustic oscillations 
(BAO) and supernovae (SN) by Komatsu et al. (2009) (black 
solid contour). The reason for this discrepancy might be that the 
ellipsoidal case describes temperature functions very well which 
are based on the mass-weighted temperature (cf. Sect. 3.2). The 
latter, however, is different from the temperature that is actually 
inferred from observations (see Sect. 3.1). A detailed study of 
this problem is therefore crucial. It will be carried out in a forth- 
coming paper. 

Besides the aforementioned contours from the potential- 
based temperature function including ellipsoidal col- 
lapse, all confidence contours are in agreement with the 
WMAP5+BAO-I-SN analysis except one: if both the temper- 
ature conversion Eq. (24) and merger effects are included in 
the calculation of the X-ray temperature function (green solid 
contour), the resulting confidence contour no longer overlaps 
with that given by Komatsu et al. (2009). 

Although we have assumed in Sect. 3.2 that the simulated 
temperature function based on Tew by Borgani et al. (2004) is 
only accidentally in good agreement with the potential-based 
temperature function including spherical collapse because of the 
overcooling problem in the central parts of simulated clusters, 
it is remarkable that the spherical function gives results that are 
compatible with results from both the WMAP5+BAO+SN anal- 
ysis and classical mass functions. At this point, it is definitely 
worth analysing whether this is only a coincidence or it really 
describes temperature functions based on observed temperatures 
better. The latter would imply that one either has to identify 
the measured temperature directly with the temperature that is 
used in the theoretical model and to include additional merger 
effects, or one has to convert the measured temperature using 
Eq. (24) to a theoretical temperature which is then used in the 
model. However, in the latter case it seems that one has to ex- 
clude merger effects since including both coiTections result in 
values for cr^ and Q,no that are inconsistent with other cosmo- 
logical probes. 

Another reason for the discrepancy between the results 
from the simulations and the comparison with cluster samples 
could be that simulated clusters are more elongated than those 
which are actually observed. The clusters for the sample by 
Vikhlinin et al. (2009a) for example are selected to appear regu- 
lar. Furthermore, the different temperature definitions are related 



in different ways to the potential shape. While Trnw traces better 
the cluster potential and hence agrees well with the temperature 
function built on ellipsoidal collapse (see Fig. 1), observed tem- 
peratures as well as Tew and as inferred from simulations 
follow the more spherical shape of the emitting gas. 

The approach including merger effects is more physically 
motivated since in the theoretical derivation of the pure X-ray 
temperature function, we explicitly assume virial equilibrium, 
which only relaxed clusters should have reached. Combining the 
results of Randall et al. (2002), who also found that mergers do 
have a significant impact on the inferred values for Q,no and erg, 
with the conclusions of Sect. 4, we believe that coiTecting for 
merger effects should be a necessary step. Note additionally that 
the scaling relation between both temperatures, Eq. (24) was es- 
tablished for clusters at z = and hence, it is not known how 
this relation evolves with redshift. 

In the preceding discussion, however, we should keep in 
mind that measurements of erg from CMB data are degenerate 
with the optical depth due to reionisation. Breaking this degen- 
eracy requires polarisation data, e.g. the T-E cross-power spec- 
trum. Thus, its value is sensitive to uncertainties in particular in 
the reionisation parameters and has changed significantly sev- 
eral times with subsequent data releases. Additional information 
from baryonic acoustic oscillations and type-la supernovae do 
not directly constrain erg either, but rather tighten constraints 
on the matter density Qmo/^^ through information on the cos- 
mological constant at fixed spatial curvature. We thus hesitate 
to accept cTg as derived from WMAP data as a firm reference. 
Weak-lensing data, that are in principle capable of constrain- 
ing cTg more directly, still yield a fairly broad range of results, 
cr% ~ 0.6 - 0.9; cf. the compilation in Bartelmann (2010). Some 
tension between expectations and data are also reflected in the 
literature. For example, while Evrard et al. (2008) prefer a high 
normalisation of the power spectrum to be consistent with nu- 
merical simulations, Reiprich (2006) concludes that data from 
the HIFLUCS sample prefer a low erg. 

Vikhlinin et al. (2009a) take merger effects into account by 
splitting the clusters of their samples into relaxed and unrelaxed 
ones by looking at their respective X-ray morphology. If a clus- 
ter is classified as unrelaxed, the mass estimate using Eq. (41) 
is multiplied by a factor of 1.17, assuming that the M-T relation 
for these two cluster samples evolves separately but similarly. 
This approach is inspired by results of a numerical simulation by 
Kravtsov et al. (2006). We think that this rigorous classification 
of clusters into relaxed and unrelaxed objects is problematic and 
should be avoided if possible. This can be done using our model 
of merger effects from Sect. 4. The resulting solid red contour 
in the lower panel of Fig. 6 is in good agreement with Vikhlinin 
et al. (2009b) (see their Fig. 3), where the rigorous classification 
was made. This and the compatibility with other cosmological 
probes indicate that our merger model can improve the determi- 
nation of cosmological parameters from X-ray data without hav- 
ing to decide individually if a cluster is relaxed or not, at least if 
mass functions and an empirical M-T relation are used to model 
an X-ray temperature function. 

7. Summary & conclusions 

In the first part of the paper, we have refined the theoretical X- 
ray temperature function of Angrick & Bartelmann (2009) in 
two different ways: First, we have used the ellipsoidal-collapse 
model by Angrick & Bartelmann (2010) to account for effects 
of the dynamics of ellipsoidal rather than spherical collapse and 
second, we have developed a simple analytic and parameter-free 
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model that takes into account the net effect of temporary X-ray 
temperature boosts of galaxy clusters that previously underwent 
mergers on the temperature function. Comparing these two mod- 
ifications to an A^-body simulation by Borgani et al. (2004), we 
have found the following results: 

- Taking into account ellipsoidal-collapse dynamics is impor- 
tant when comparing the theoretical model to temperature 
functions of numerical simulations that are based on mass- 
weighted temperatures averaged over a large volume (e.g. in- 
side the virial radius). This is the temperature definition con- 
sistent with the virial theorem, Eq. (11). Temperatures from 
real observations, however, are similar to the spectroscopic- 
like temperature Tsi or the emission-weighted temperature 
Jew We have shown that temperature functions based on 
Jew are in good agreement with the theoretical model only 
for spherical-collapse dynamics. We have guessed that this 
is mainly due to the overcooling problem in the centres of 
simulated clusters and shifts the measured cluster temper- 
atures to larger values. Therefore, it is most prominent for 
temperatures inferred from the central cluster regions. 

- Especially for z > 0.5, the effects of mergers cannot be 
excluded since the higher the redshift, the more clusters 
are unrelaxed and therefore deviate from virial equilibrium. 
Our simple analytic and parameter-free model based on the 
merger probability derived by Lacey & Cole (1993) can ac- 
count for these effects. Including it in our theoretical mod- 
elling leads to a substantially improved agreement with the 
simulation for redshifts < z < 1 . 

In the second part, we have used both mass-based and 
potential-based X-ray temperature functions including either 
spherical- or ellipsoidal-collapse dynamics together with sam- 
ples of Ikebe et al. (2002) and Vikhlinin et al. (2009a) to con- 
strain the cosmological parameters Qmo and erg. We have anal- 
ysed the influence of merger effects on the inferred values of 
both parameters by using our analytical model. In addition, we 
have tested whether it is necessary to convert the measured tem- 
peratures by means of Eq. (24) or if it is possible to compare 
them to the theoretical prediction directly. The main results are 
the following: 

- The temperature function based on spherical-collapse dy- 
namics by Angrick & Bartelmann (2009) leads to confidence 
contours in the Q,„o-C8 plane that are compatible to those in- 
ferred from classical mass-based temperature functions. The 
best-fit values might be shifted by a few percent. However, 
the formalism does not refer to cluster masses whose relation 
to observables have to be calibrated. 

- Using the potential-based temperature function that incorpo- 
rates ellipsoidal collapse as the theoretical model yields con- 
fidence contours that are not compatible with a joint analy- 
sis by Komatsu et al. (2009) using WMAP5+BAO+SN data 
since the inferred values for Q.^o seem too large. We believe 
that the reason is the discrepancy between mass-weighted 
temperature to be used in the virial theorem and the actual 
measured temperature that is similar to Tsj, but also selec- 
tion effects such as preferring more regular clusters over el- 
lipsoidal clusters in a cluster sample as well as different re- 
lations of the various temperatures to the potential structure 
(Jmw traces more the cluster potential, whereas Tgw and T^i 
follow mainly the spherical shape of the emitting gas) in- 
troduce a bias that can lead to the discrepancy between our 
analysis based on the temperature function incorporating el- 
lipsoidal collapse and the WMAP5-I-BAO-I-SN data. 



- In case of the temperature function built on spherical col- 
lapse, identifying the measured temperature directly with the 
temperature that is used in the theoretical model and ad- 
ditionally including merger effects, or converting the mea- 
sured temperature using Eq. (24) before comparing to the 
model and disregarding merger effects give similar results 
for Qn,o and erg. Since the comparison to the simulation and 
the results from the mass-based temperature functions indi- 
cate that merger effects lead to biases, we believe that the di- 
rect comparison between measured and theoretical temper- 
atures is the correct choice when using our potential-based 
X-ray temperature function. 

- Almost all of our results using the potential-based temper- 
ature function relying on spherical collapse are compatible 
with the WMAP5+BAO+SN data. Only if both the temper- 
ature conversion Eq. (24) and merger effects are taken into 
account, the resulting confidence contours disagree with the 
latter. Since merger effects do have a significant impact on 
the determination of Qn,o and erg (one of the main results 
of Randall et al. 2002), a temperature conversion does not 
seem to be necessary in the context of the potential-based 
temperature function. The reason why the spherical temper- 
ature function gives such good results at all has to be further 
analysed since the agreement with the data from the simula- 
tion of Borgani et al. (2004) is suspected to be mainly due to 
the overcooling problem in the centre of simulated clusters. 
A deeper understanding of how to relate measured tempera- 
tures to those used in our theoretical models is necessary. 

- The combination of the mass function by Tinker et al. (2008) 
with our simple analytic merger model yields similar results 
as the technique by Vikhlinin et al. (2009a,b), i.e. rigorously 
classifying galaxy clusters into relaxed and unrelaxed ob- 
jects. However, we believe that our model is physically better 
justified since we take merger effects into account statisti- 
cally. 

- Although the results for the samples by Ikebe et al. (2002) 
and Vikhlinin et al. (2009a) are in agreement with each other, 
pronounced differences exist, implying that the final confi- 
dence regions from our analysis of the data might still be 
biased from systematics in the X-ray analysis and therefore, 
even better agreement with other cosmological probes could 
be achieved. 
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